b = 0.02;
I = 0.04;
k = -1000;

roots([1 b/I 0 k/I])
